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Theoretical  Analysis  of  Image  Processing  Using  Parameter-Tuning 

Stochastic  Resonance  Technique 

Bohou  Xu,  Zhong-Ping  Jiang,  Xingxing  Wu,  and  Daniel  W.  Repperger 


Abstract — Parameter-tuning  stochastic  resonance  has  been 
successfully  applied  to  the  one-dimensional  signal  processing. 
This  paper  explores  the  feasibility  to  extend  this  technique 
for  image  processing.  Based  on  the  two-dimensional  nonlinear 
bistable  dynamic  system,  the  equation  satisfied  by  the  system 
output  probability  density  function  is  derived  for  the  first  time. 
The  corresponding  equation  for  the  one-dimensional  system 
is  the  famous  Fokker-Planck-Kolmogorov  (FPK)  equation.  The 
stationary  solution,  eigenvalues  and  eigenfunctions  of  this  equa¬ 
tion  are  then  investigated.  The  upper  bound  of  the  system 
response  speed  and  the  related  calculation  algorithm  which 
are  necessary  for  the  applications  of  this  technique  to  image 
processing  are  also  proposed  in  this  paper.  Finally,  the  potential 
applications  of  this  approach  in  image  processing  and  some 
future  research  are  suggested. 

Index  Terms — Stochastic  Systems,  Stochastic  Resonance,  Fil¬ 
tering,  Nonlinear  Systems,  Image  Processing 

I.  INTRODUCTION 

Image  processing  has  been  widely  applied  in  different 
areas,  such  as  diagnosing  tumors  in  medical  images,  de¬ 
tecting  and  identifying  hostile  targets  in  military  images. 
Over  the  years,  many  effective  image  processing  algorithms 
have  been  proposed  to  meet  the  increasing  requirements  on 
image  qualities.  For  the  images  corrupted  by  noise,  most  of 
the  denosing  algorithms  will  try  to  remove  or  suppress  the 
noise  from  the  systems,  because  the  noise  is  usually  thought 
to  be  annoying.  Stochastic  resonance,  on  the  contrary,  is  a 
phenomenon  that  the  noise  can  be  used  to  enhance  rather  than 
hinder  the  system  performance.  The  concept  of  stochastic 
resonance  was  first  proposed  by  Bcnzi  in  1981  [11.  Since 
then,  stochastic  resonance  has  been  applied  in  a  wide-range 
of  areas,  such  as  physics,  chemistry,  biomedical  sciences,  and 
engineering  [2],  It  has  been  successfully  used  to  improve 
the  balance  control  for  elderly  people  [3].  The  profoundly 
deaf  people  can  improve  their  speech  understanding  with 
the  aid  of  noise  [4].  For  the  signal  processing  area,'  the 
stochastic  resonance  technique  has  been  applied  to  the  signal 
detection  [5],  signal  transmission  [6],  and  signal  estimation 
[7].  In  order  to  realize  stochastic  resonance  to  make  the  noise 
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beneficial  to  the  system,  the  synchronization  between  the 
input  signal  and  the  noise  must  occur.  Basically,  there  are 
two  approaches  to  realize  the  stochastic  resonance  effect. 
The  traditional  method  is  to  add  an  optimal  amount  of  noise 
into  the  system  [2].  Parameter-tuning  stochastic  resonance 
proposed  by  us  is  the  other  way  f 8 ]-[  1 1  ][  1 6].  It  realizes 
the  stochastic  resonance  effect  by  timing  system  parameters 
to  their  optimal  values  without  adding  any  additional  noise 
into  the  system.  We  also  reveal  that  the  parameter- tuning 
stochastic  resonance  is  superior  to  the  traditional  method  [8], 
especially  when  the  initial  noise  intensity  is  already  beyond 
its  resonant  region.  Parameter-tuning  stochastic  resonance 
has  also  been  used  to  recover  the  noisy  multi-frequency  sig¬ 
nals  [9],  reduce  the  bit-error  rate  (BER)  of  the  transmission 
of  baseband  binary  signals  [8].  Image  processing  is  another 
potential  application  area  of  stochastic  resonance  technique. 
There  are  some  initial  research  work  on  it  [I2][I3].  The 
initial  research  results  demonstrate  that  it  is  promising 
and  feasible  to  develop  the  innovative  and  effective  image 
processing  algorithms  using  stochastic  resonance  technique. 
All  these  approaches,  however,  are  based  on  simulations 
and  are  lacking  rigorous  theoretical  analysis.  It  is  difficult 
to  implement  and  extend  to  other  image  processing  fields. 
Another  problem  with  these  approaches  is  that  they  need 
to  add  an  optimal  additional  noise  into  the  images,  which 
is  impossible  for  some  image  processing  tasks.  All  these 
motivate  us  to  investigate  the  application  of  parameter- 
tuning  stochastic  resonance  in  image  processing  based  on  the 
systematic  and  theoretical  analysis.  This  paper  will  mainly 
focus  on  the  theoretical  investigation  of  the  feasibility  of 
this  approach.  In  order  to  apply  the  current  parameter- 
tuning  technique,  the  two-dimensional  image  signals  can 
be  first  converted  to  onc-dimcnsional  signals.  Unfortunately, 
our  research  demonstrates  this  method  is  not  effective.  The 
only  possible  method  is  then  to  treat  the  two-dimensional 
image  signals  directly  without  conversion  and  develop  all 
the  theoretical  results  in  an  analytic  way  as  one-dimensional 
parameter-tuning  stochastic  resonance.  For  one-dimensional 
systems,  the  parameter-tuning  stochastic  resonance  is  based 
on  the  derivation  of  the  Fokker-Planck-Kolmogorov  (FPK) 
equation  satisfied  by  the  system  output  probability  density 
function  and  also  the  derivation  of  the  system  response  speed 
and  its  calculation  algorithm.  Similarly,  for  two-dimensional 
image  signals,  the  feasibility  investigation  means  whether 
the  equation  and  its  solutions  satisfied  by  system  output 
probability  density  function  can  be  derived  and  whether  the 
system  response  time  can  be  calculated.  So  far,  no  related 
research  result  could  be  found  in  this  area, 


The  rest  of  this  paper  is  organized  as  follows.  In  Section 
IT,  we  will  first  propose  a  two-dimensional  nonlinear  bistable 
dynamic  system  and  then  derive  the  equation  satisfied  by 
the  system  output  probability  density  function.  Section  III 
will  give  the  stationary  solution  of  this  equation.  The  system 
response  speed  for  this  two-dimensional  system  will  be 
investigated  in  Section  IV.  The  potential  applications  of  this 
approach  in  image  processing  are  mentioned  in  Section  V. 
Finally,  Section  VI  closes  the  paper  with  brief  concluding 
remarks  and  future  research  directions. 

II.  TWO-DIMENSIONAL  NONLINEAR  BISTABLE 
SYSTEM  AND  RELATED  FPK  EQUATION 


where 

M0=  1,  (6) 

Mn{u,  x,  Ax,  y,  Ay) 

(v  —  u)n  P(v,  x  +  Ax.  y  +  Ay\u.  x .  y)dv ,  n  =1.2,  ...(7) 

Then 

p  (w,  x  +  Ax,  y  +  Ay)  —  p(w,  x,  y) 

=  ^>  — j  /  M-n  (u,  x,  Ax,  y,  Ay)p(u,  x,  y) 

n= 1  ’  ‘ 

re  1  1 

=  £  ~(--^)n[Mn{iu,  x.  Ax,  y,  A y)p{w,  x.  y)}.  (8) 

n  —  i 


The  one-dimensional  nonlinear  bistable  stochastic  reso¬ 
nance  system  can  be  described  by  the  following  equation 
[2] 

x(t)  =  ax(t)  -  bx3(t)  +  .<?(£)  +  rj(t),  (1) 


where  a  and  b  are  system  parameters,  y(t)  is  the  input  signal, 
and  7](t)  is  an  additive  Gaussian  white  noise  with  zero  mean 
average  and  autocorrelation  of  {//(£) 7/(0))  =  2 DS(t). 

For  this  nonlinear  dynamic  system,  the  output  signal-to- 
noise  ratio  will  be  maximized  when  an  optimal  amount 
of  noise  is  added  into  the  system.  This  is  the  stochastic 
resonance  phenomenon. 

Similarly,  we  can  propose  the  following  two-dimensional 
nonlinear  bistable  dynamic  system 


d2tv 

dxdy 


f(w)  +  F  ( x  ,y). 


(2) 


where  w  =  w{x,y)  is  the  state  variable  (system  output), 
f(w)  =  aw  —  bw3  +  h,  r (x,y)  is  additive  white  Gaussian 
noise,  and  h  is  the  input  signal. 

The  corresponding  difference  equation  is 


w(x  +  Ax,  y  +  Ay)  ~w(x  +  Ax,  y)  +  vj(x,  y  +  Ay)  —  w(x,  y) 
+,/'('w)A.t:A  y  +  r(x,  y)  Ax  Ay,  (3) 


In  order  to  derive  the  FPK  equation  satisfied  by  the  system 
output  probability  density  function,  one  possible  approach  is 
to  use  the  similar  method  as  deriving  one-dimensional  FPK 
equation  [141. 

We  have 


p(w ,  x  +  Ax,  y  +  Ay) 


w ,  x  +  Ax,  y  +  Ay\u,  x,  y)p(u.,  x ,  y)du, 


(4) 


where  p{w,  x,  y)  is  the  probability  density  function  of  the 
system  output  w  and  P{w,x  +  A x,y  +  Ay\u,x,y)  is  the 
conditional  probability  density  function. 

Also, 


P(w,  x  -f  Ax,  y  +  Ay  [a,  x ,  y) 

=  j  S ( v  —  w ) P ( v ,  x  -f  A x ,  y  +  A y (it, M  y ) d\ 
(v  -u)n  dnS(u-iv) 


=  E 


dip- 


P(v,  x  +  A  x .  y  -f  A  y  |  u ,  x ,  y )  dv 


Similarly,  we  can  obtain 

p(w  ,  x  +  Ax,  y)  —  p(w ,  x.  y) 

oc 

=  L  —(“Tr  X.  Ax.  y.  0  )p(w.  ./■.  ;</)!•  (9) 


p(;w  ,  x,  y  +  Ay)  -  p(w.  x,  y) 

i(-A 

n! v  dw 


00  i  d 

=  y,  X,  0.  y,  A y)p(w,  x,  yy)].  (10) 


From  this,  we  can  derive  the  corresponding  FPK  equation 

d2p(w.x.y)  1  .  d  w,  r„  .  .  . 

=  [C„(w,x,w)p(w,a:,y)),(ll) 


dxdy 


where 

C„(u),x.y)=  1  i in  —  [ Mn  ( w .  x,  Ax,  y .  Ay) 

' a.t-.o  Ax  Ay  v 

Ay->0 

—  M a  (w ,  x ,  A x ,  y ,  0 )  —  Mn  (w ,  x ,  0 ,  y ,  Ay )  ] .  ( 1 2) 


For  n  =  1,  C\  can  be  calculated  as 

Ci  (u>,  x.  y)AxAy=(vj(x  +  Ax.  y  +  Ay)  -  w(x  +  Ax,  y) 
-w(x,  y  +  Ay)  +  w{x.  y)) 

=f(w)AxAy,  ( 1  3) 

where  the  angular  brackets  denote  Lhe  ensemble  average. 

Unfortunately,  it  is  impossible  to  calculate  C„  for  n  >  2, 
because 

Cd(w,  x,y)AxAy—([w(x  X  A:r.  y  t  A//)  --  u  ( x.  y/)]2) 

—  { [ta(x  -I-  Ax.  y)  -  w{x.  y) j“) 

-  ( [x;  ( x .  y  H-  A  y )  -  it 1  ( x .  .{/)  J J ) .  (14) 


The  system  described  by  (2)  provides  no  enough  informa¬ 
tion  to  calculate  the  Cn  for  n  >  2. 

To  overcome  the  above  difficulty,  we  reduce  the  partial 
differential  equation  (2)  to  an  ordinary  differential  equation 
along  the  line  x  =  x0  -f  tAx,  y  =  y0  -f  t  Ay  as 

d2w 

=  AxAyj(w)  -f-  AxAyTixo  +  (Ax,  y(]  +  tAy),  (15) 

where 


oo 


(r(x,7/)r(.r ,.?/,))  =  2 DS(x  -  :r,.y-  yi). 


(16) 


The  equaLion  (15)  can  then  be  rewritten  as 


dw  . 

~—=Axv, 

dt 

^=Ayf(w)  +  AyT(Axt,  A yt) 
=Ayf(w)  +  AyTi(t). 


Assume  po(w,v)  =  e  a°v  we  have 

—Axve~a'°0  ip  (to)  +  'yAye~a°v  y>{w) 

+ 2 au A y ( -  •  7W-b  f (w)) v e~a°v  ip{w) 

+DAy2(  —  2aQ  +  4o,QV2)e~a{)V  y>(w)  =  0.  (23) 


where 


-2 D8(t-ti).  (18) 


Based  on  the  concept  in  [14],  wc  can  prove  the  proba¬ 
bility  density  function  p(w,v,t)  satisfies  the  following  FPK 
equation 

dpiw^t)  =- A[Ax-tjp(-u; ,v ,/.)]  -  2L[Ayf{.vj)p(w,v,t)\ 

+  (19) 
avz 

Equation  (19)  is  independent  on  (x0,;(/o)  explicitly,  the 
solution  of  this  equation  will  be  valid  for  any  (xo,yo)- 

Equation  (19),  however,  docs  not  have  a  stationary  solu¬ 
tion,  because  it  is  a  hyperbolic  equation  without  any  damping 
[17].  In  order  to  overcome  this  problem,  equation  (2)  is  now 
changed  to 


In  order  to  meet  (23),  the  following  should  be  satisfied  by 
the  coefficient  of  v\  for  i  —  0.  1, 2 

a0:  7  A  y  —  2aoDAy2  =  0.  (24) 

a1:  —Axy:  (  w)  -(-  2o.oAyf(tu)yj(w)  =  0.  (25) 

v2:  —2  any  Ay  -f-  Aa'^DAy2  —  0.  (26) 


By  solving  (24),  (25),  and  (26),  we  can  get 


aa  = 

7 

2D  Ay'’ 

(27) 

<p{w) 

=  /Voexp[^ 

-  f'  /('«: 

Jo 

OH 

(28) 

po(w, 

v)  =  No  exp[- 

7  V 2 

2D  Ay 

+  DAxJ0  1 

[w)dw],  (29) 

•e  N0 

can  be  determined  by 

r  ~r  O «  + 

J  _  cc  J  —  r 

■  OO 

pu{vx  i 

>)dwdv  ~  1. 

(30) 

d2-w 

dxdy 


+  f(w)  +  r(sc,y), 


(20) 


where  7  is  a  positive  damping  coefficient. 

In  the  similar  way,  we  can  derive  the  FPK  equation  for 
system  (20) 


dp(w,  v.  t) 
dt 


=-A[^XVp(u,,v.t)\ 

-^l&yi-rv  +  f(v>))p(w,v,  t)) 


+  DAy 


■2  d2p(w,v,t) 
dv2 


(21) 


In  the  following  section,  we  will  show  that  the  stationary 
solution  of  (21)  exists.  System  described  by  (20)  will  be 
used  as  the  two-dimensional  bistable  dynamic  stochastic 
resonance  system  to  process  the  image  signals. 


For  convenience,  (29)  can  be  rewritten  as 

po(w.v)  = 

where 


(31 


■  /  >  7  2  7 

0(w.v)  =  — v 


2D  Ay  DAx  J0 


f(w)dw  -  In  .V().  (32) 


IV.  SYSTEM  RESPONSE  SPEED  OF 
TWO-DIMENSIONAL  BISTABLE  SYSTEMS 

Similar  to  one-dimensional  parameter-tuning  stochastic 
resonance,  investigating  the  characteristics  of  system  re¬ 
sponse  speed  and  developing  its  calculation  algorithm  are 
very  important  tasks  for  the  two-dimensional  parameter¬ 
tuning  stochastic  resonance.  It  will  determine  whether  the 
one-dimensional  parameter-tuning  stochastic  resonance  can 
be  extended  to  the  two-dimensional  case. 

Let 


III,  STATIONARY  SOLUTION  OF 
TWO-DIMENSIONAL  FPK  EQUATIONS 

Equation  (21)  is  linear  for  the  probability  density  function 
p(w,v,t),  so  it  is  possible  to  be  solved  by  eigenfunction 
expanding  method. 

Let  po(w,v)  be  the  stationary  solution.  It  will  satisfy  the 
following  equation 

-  [Axvpo (w,  t’)]“  [ Ay ( -yv  +  f(w))p0(w,  v)] 

+m,4>(7")=0,  (22) 

dv2 


Dy—V  Ax, 

D'2~Ay[~yv  +  f(w) 
D2-2  =  DAy2. 

Equation  (21)  can  be  rewritten  as 
d 


dt  Oil 
=Lpr(p) 


,  n  ,  d  ,  „  .  „  02(> 

(■ Dip )  ~  a2,(°2p)  + 1)12 oa 


where 


<P)  = 


(33) 

(34) 

(35) 


(36) 


(22) 


(37) 


Assume  p  =  £(w}v)e  xt,  equation  (36)  then  becomes 

A  £(w,v)  =  -Lfp£(w,v).  (38) 


Now,  we  can  calculate  the  symmetric  part  Ls  and  the  anti¬ 
symmetric  part  Because 


Similar  to  the  one-dimensional  case,  let 

£{w,v)  —  'ip(w,  v)e~(fl/2, 


(39) 


where  0  is  defined  in  (32)  and  e  ^  is  the  stationary  solution 
of  (2 IX 

In  this  case,  equation  (38)  becomes 


Xxjj  =  —  7/0. 

where  L  is  a  differential  operator,  and 

L'ip=e^2  Lfp(£-^/2,4>) 


(40) 


ov 


dwK 
+eiinD2 

ovz 


(41) 


We  can  assume 

f(w,v)  =  0(e_<^2),  as  w  — »  ±oo  or  v  — »  ±oo.  (42) 
because 


-  lL(D2e:  *)  +  ±(Dn$-<T*\  =  0.  (49) 
ow  ov  ov  Ov 


we  can  get 


r  -  (L+LC) 
LsV=- - o - 


1 , 


J>j'l  0  f  n  -0/2^  J>/2  & 


+e*nDi2^(e~*n1>)  +  <TmDx  f(e03ib) 

OV 2  OV) 


-he 


lD2~{e0/1V>)+e- 
ov 


r  <r  -  . 
Iho—{*  ) 

(Jv~ 


04>n  d 


(50) 


lim  p(w,v,t.)  =  lim  p(w.vpt)  =  0,  (43) 

'  — >  ±00  v— *±cx; 


e  ^  diodv=const  0. 


(44) 


Also 


We  denote  the  conjugate  operator  of  L  as  Lc .  Let  Ls  = 
(. L  +  Lc)/ 2,  and  Las  =  (L  —  Lc)/2 ,  that  is  L  =  Ls  +  Lns, 
we  can  obtain  the  following  theorems. 

Theorem  I:  The  differential  operator  L  in  (40)  can  be 
decomposed  into  symmetric  part  L#  and  anti-symmetric  part 
Las.  Also,  Ls  is  negative  semi-definite. 

Proof:  First,  we  will  calculate  the  conjugate  operator 
of  L.  It  is  denoted  as  Lc.  Assume 


La»i/>=-(L  -  L(:)ij>  —  (L  —  L,)v 

=e*/2  {-f(Die~,unii>\-5L(p2e--  ",/ 2 ,,, ) 


d_ 

6  dv K 


=  -e*/2{-f-(Die-*/2iJ>)  +  ™(D2e"(/V2V;) 

Ow  ov 

i  ^  rn  ^0  —  0/2  ji] 

+  lTv{D22Tve  V;)}- 


(51) 


lim  +  (+;,  v)  =  lim  rf  tu.v)  =  0,  (45) 

t/j— *±oo  w  —  *±oc 


lim  '0  ( w ,  v )  =  lim  q  (w ,  i? )  =  0 . 

*±oo  tj— >±oo 


(46) 


From 

rj  Ldv  tpdw 


//- 


,  0/2  n  ^  /  —0/2 


■0 )  ]  r7«;  dv 


-  j  j  0[^i  e  ' f'V2  ^  (e''V2?/)  +  L»2c  4!/2^vf/!j/'2v) 

+  £>22  2  ( c r>/  2  r/ )  ]  r/iedt; , 

we  can  define  the  conjugate  operator  of  L  as 


(47) 


(48) 


To  conclude,  we  get 


0/2  ^ 


Ls=e'p/~ \D-22C  ~'l> d-(e''^~)\. 

Ov  Ov 

La,=-e*r2\d-{nie-<!'/2)  +  f  {D2e  */2) 

OW  Ov 

,  .//)  ^0  0/2 ^ 


'dv1  dv 1 

=  j  j  e,l,/2^f[D22e~':‘ ~  (et,/2n)\dwdv 
=  /  j  4!£-sr/+w.’d'?;.  (54) 


(52) 


(53 1 


For  any  0  and  7j  satisfying  (45)  and  (46),  wc  have 
j  j  T)Lsij)dxudv 

=  lf  e*'2vllDne-*£.y'**)\du>dv 


and 


I  j  r}Las'4}dwdv 


+  -^-(D22^e  'i,/2ip]\d;wdv 


=  //  le  ^/2^Di  (e0/2r/) +  e‘  */2^2 1; (e*/2/?) 

+D22  T^-G~q>l/2'fpS-  ((d^2rj)\du)dv 


'  Ac 


Ac 


±~{D22^-e  <b/2i t)]dwdv 
dv  dv 


=  ~  j' J  d}Las't 


77  dwdv. 


(55) 


This  means  the  operator  Ls  is  symmetric  and  Las  is  anti¬ 
symmetric.  Now,  we  will  prove  Ls  is  also  negative  semi- 
definite. 

We  have 

j  j  ipL s  f dw dv= j  j  ip e  ■' /2  ~  [ £>22 e ^  ^  ^  ( e 2 01 ) ]  dw dv 

D22&~^ {-^~^^2dA2 diudv  <  0.  (56) 
ov 


From  this,  we  can  conclude  that  Ls  is  semi-negative 
definite.  This  completes  the  proof.  ■ 


Substitute  (61)  into  (60),  we  get 

,  If  E^o  EiEo  Cji^i  LjilCjkdwdv 

YZo\Cn\2  '  (  ) 

We  can  also  obtain  the  following 

EE  Cj'Cjk  f  I  d>*  Li>id:wd.v 

:-0  k=Q  ’  ■’ 

=  X]  ^  ^(~'Kr^rk+  j  j  Ln  -  ?/  0’  (ft*  v/f ’ ) .  ( 63  ) 


y  ,  A j.CjkMd  LasTpk  +  /A  CjjCjkP J  Ltaaty'k 
/-0  Jt  =  0  1=0  k=0 

OC  OO  OO  DC 

=EE  CjiCjk'lpt  Tasfk  T  ^  ^  ^  ^  CjiCjkdd  Tafifk 


7-^0  A*=0 
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From  above,  we  have 


ncnt)  TZo\C^ 

(  j)  YZ«\Cp? 

Assume  0  <  /tle(Ai)  <  /?e(A2)  <  ...,  then 

Rc(x  )  =  S”|  MX  ^ 

Rc(\y)  ^A" 

Here,  Cw  =  0.  We  then  prove  that  Ref A)  >  A*. 


(64) 


(65) 


(66) 


Theorem  2:  If  A  is  any  non -zero  eigenvalue  of  the  op¬ 
erator  —  L,  and  A f  is  the  least  non-zero  eigenvalue  of  the 
operator  -L, 3,  then  Re  A  >  Af. 

Proof:  Assume  { A,f ,  ft }  ,  for  i  —  0,1,...,  are  the 
eigenvalues  and  eigenfunctions  of  —  L*,  that  is 

-L^fl  =  Xfff  i  =  0,  1,  ...  (57) 


Because  —  Ls  is  symmetric  and  semi-positive  definite,  we 
have 


0  =  AS  <  A,  <  A5  < 


iptAjd/wdu 


A,. 


(58) 

(59) 


Moreover,  we  assume  that  Ay  and  fj  are  the  eigenvalue 
and  eigenfunction  of  —L.  That  is  —  Lfj  =  A jfj.  \3  might 
be  complex.  Then,  we  get 


fj  1 pj  Lip  j  dwdv 
If  fj  fj  dwdv 


(60) 


where  -f*  is  the  complex  conjugate  of  f  j. 

Because  {t/A}  is  complete  in  H-space,  fj  can  be  expanded 
as 


*  =  E^'  (61) 

7=0 


From  Theorem  2,  we  can  regard  Aj  as  a  lower  bound 
of  the  system  response  speed  of  —  L  and  take  X\  as  an 
approximation  of  the  system  response  speed  of  system  (20). 
We  now  investigate  its  calculation  algorithm. 

The  eigenvalue  problem  of  operator  Ls,  AN/;  =  -  LsiP.  can 
be  rewritten  into  the  variational  form  (Rayleigh  quotient)! !  5] 


A’,=  st. 


=  St. 

v^o 


=  St. 

V.^0 


=  st. 

C4n 


JJ  fLsijxlwdv 
If  f  2  dwdv 

ff  ip2  dwdv 

n_D2 

f  f  f 2 dwdv 

.jJjV7e^(§fdwdv 
ff  e-^Cfdwdv 


(67) 


where  0  =  e^2'?/;,  and  “st.”  denotes  the  stationary  value  of 
the  functional. 

Obviously,  As  >  0.  It  will  take  the  zero  value,  if  C  is 
constant  (This  means  f  =  Ce-^2,  and  it  is  a  stationary  so¬ 
lution).  The  above  problem  can  be  solved  using  the  subspace 
iteration  method. 

The  calculation  algorithm  of  Aj  is  outlined  as  follows 

1)  Determine  the  integral  range  in  (67); 

2)  Divide  the  integral  range  into  n  x  m  parts  and  use 
the  finite  element  method  to  transform  the  variational 
problem  (67)  into  a  general-matrix  eigenvalue  problem 


(\[M}-[K}){a  =0. 


(68) 


where  [M]  is  a  positive  definite  matrix  and  [K]  is  a 
non-negative  matrix. 

3)  Use  subspace  against-iteration  method  to  get  the  min¬ 
imum  positive  eigenvalue  A*. 

Besides  the  stationary  solution,  we  can  now  get  the 
expanding  solution  of  (21) 

p(w,V .  t)  =  po(w,v)  +  Cl1pi(w,v)c~Xlt'-0^ 

+C'2'02 (w ,  v)e ~ X2f  ~ 0/2  +  . . .  (69) 

where  p$  is  the  stationary  solution  described  by  (31),  and 
{A for  i  =  1,2,...,  are  the  eigenvalues  and  eigen¬ 
function  of  —L.  It  is  obvious  that  p(w,v,t)  &  p$(w,v), 
when  A \t  »  1.  This  can  be  satisfied  by  adjusting  the  system 
parameters  a  and  b  which  are  defined  in  (2).  In  this  case,  we 
can  use  pu(wpu)  to  approximate  p(w,v,t)  when  applying 
this  technique  to  image  processing. 

V.  POTENTIAL  APPLICATIONS  IN  IMAGE 
PROCESSING 

The  derivation  of  (69)  demonstrates  the  feasibility  to 
extend  the  concepts  of  parameter-tuning  stochastic  resonance 
to  the  two-dimensional  case.  The  two-dimensional  bistable 
system  (20)  can  be  used  as  a  nonlinear  filter  to  process  the 
noisy  two-dimensional  image  signal  to  improve  the  image 
quality.  The  probability  characteristics  of  the  image  signal 
after  the  filter  is  described  by  (20).  It  can  then  be  used 
to  derive  other  performance  measures  based  on  different 
image  processing  tasks,  such  as  image  signal-to-noisc  ratio, 
probability  of  target  detection  error,  etc.  The  performance 
measures  can  then  be  taken  as  the  object  function  to  be 
optimized  by  tuning  the  system  parameter  a  and  b,  which  is 
the  concept  of  parameter-tuning  stochastic  resonance.  This 
method  has  the  potential  to  be  applied  in  image  processing 
tasks  in  which  the  image  signals  are  corrupted  by  the  noise. 
The  system  parameters  will  be  tuned  properly  to  synchronize 
the  signals  and  the  noise  to  convert  the  noise  to  be  a  positive 
factor  to  improve  the  image  qualities. 

VI.  CONCLUSION  AND  FUTURE  WORK 

Through  theoretical  analysis,  this  paper  reveals  that  the 
one-dimensional  parameter-tuning  stochastic  resonance  can 
be  extended  to  the  two-dimensional  case  and  it  is  feasible  to 
use  it  for  image  processing.  The  Fokker- PI  a  nek -Kolmogorov 
(FPK)  equation  for  the  two-dimensional  nonlinear  bistable 
dynamic  stochastic  resonance  system  can  be  derived.  Its 
stationary  solution  is  also  given  in  this  paper.  Also,  the  sys¬ 
tem  response  speed  for  this  two-dimensional  bistable  system, 
which  is  another  important  part  of  parameter- tuning  stochas¬ 
tic  resonance,  can  be  calculated  with  our  proposed  algorithm. 
All  these  demonstrate  that  all  the  critical  techniques  of 
the  parameter-tuning  stochastic  resonance  are  now  available. 
They  are  ready  to  be  applied  to  process  two-dimensional  im¬ 
age  signals.  The  two-dimensional  parameter-tuning  stochas¬ 
tic  resonance  provides  an  innovative  and  promising  approach 


for  image  processing.  It  will  find  wide-range  applications, 
just  like  that  of  one-dimensional  parameter-tuning  stochastic 
resonance.  Next,  we  will  apply  this  new  approach  to  real 
image  processing  tasks.  We  will  start  from  the  binary  images 
with  white  Gaussian  noise  and  then  extend  to  gray  images 
with  colored  noise. 
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